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We report on the magnetic trapping of an ultracold ensemble of ®^Rb atoms close to a super¬ 
conducting ring prepared in different states of qnantized magnetic flux. The niobinm ring of 10 pm 
radins is prepared in a flux state n<l>o, with $o = h/2e the flux quantum and n varies between ±5. 
An atomic cloud of 250 nK temperature is positioned with a harmonic magnetic trapping potential 
at ~18pm distance below the ring. The inhomogeneous magnetic field of the supercurrent in the 
ring contributes to the magnetic trapping potential of the cloud. The induced deformation of the 
magnetic trap impacts the shape of the cloud, the number of trapped atoms as well as the center- 
of-mass oscillation frequency of Bose-Einstein condensates. When the held applied during cooldown 
of the chip is varied, the change of these properties shows discrete steps that quantitatively match 
flux quantization. 

PACS numbers: 37.10.Gh, 74.25.Ha 


The coherent coupling between atoms and single flux 
quanta in a superconducting circuit is an important 
ingredient of future cold atom-superconductor hybrid 
quantum systems in which quantum states are trans¬ 
ferred from one system to the other. The construction of 
such a hybrid quantum system is targeted in a number 
of recent experiments and proposals [IHS], and should 
allow the study of fundamental interactions between the 
two systems pTl4T^ . One of the most prominent conse¬ 
quences of the existence of a macroscopic wave function 
in superconductors is the quantization of the magnetic 
fluxoid, which has been shown e.g. in superconducting 
rings and cylinders Eoi m]. In superconducting atom 
chip experiments, trapped Abrikosov vortices have been 
used to magnetically trap atoms in spatially inhomoge¬ 
neous fields P^27) . These traps are affected by the mo¬ 
tion of the vortices that potentially cause heating and 
losses of the cold atoms |5S1 [H]. Pinning the vortices 
would suppress this noise source and could be used to 
generate subwavelength magnetic lattices [30] , as well as 
hybrid quantum systems based on atom traps formed by 
single pinned flux quanta m- The creation of a flux 
superposition state could give rise to a superposition of 
the magnetic trapping potential and therefore of the po¬ 
sition of an atomic ensemble [32]. Ultrafast coupling (10 
ns) between cold Rydberg atoms and SQUIDs has also 
been theoretically predicted [35]. It is therefore crucial 
to understand the impact of single flux quanta onto an 
ensemble of trapped ultracold atoms. 

In this Letter, we report on how a discrete number of 
flux quanta stored in a superconducting ring affects the 
trapping parameters of a superimposed magnetic trap. 
The discrete nature of the magnetic flux in the ring is 
observed both in the atom number and the oscillation 



FIG. 1. (a) Sketch of an atomic ensemble trapped at the su¬ 
perconducting ring, (b) Optical image of parts of the super¬ 
conducting chip. Shown are four trapping wires, the 100 pm 
wide trapping wire used for the present experiment is high¬ 
lighted in blue. The location of the confinement wire is 
sketched in yellow. The center of the superconducting ring 
is located 70 pm from the right edge of the trapping wire. 
The chip is mounted upside down on the cryostat, so gravity 
points in +z direction. 


frequency of atoms in the trap. 

We magnetically trap an ensemble of cold ®^Rb atoms 
on a superconducting atom chip and guide it to the 
vicinity of a superconducting ring, as sketched in Fig. 
[^a). The atom chip (Fig. [^b)) is a sapphire substrate 
(~330pm thick) with patterned niobium thin film struc¬ 
tures (thickness d = 500 nm). The chip contains several 
Z—shaped lines (“trapping wires”) of different widths 
used for trapping and moving the atoms, the broadest 
of which (100 pm wide) is used to trap atoms in the ex¬ 
periments described in this Letter. The ring has an in¬ 
ner radius n =9 pm and an outer radius Tq =11 pm. It 
is placed 70 pm from the edge of the trapping wire in y 
direction. The superconducting atom chip [^ is attached 
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to the cold finger of a helium flow cryostat at tempera¬ 
ture T = 4.2 K. The atoms are prepared in the hyperfine 
ground state 5 S 1 / 2 F = 1,tof = —1 in a room temper¬ 
ature part of the setup and subsequently transported to 
a position below the superconducting chip by means of 
optical tweezers, as detailed in Ref. |34j . 

The microtrap is realized by the superposition of the 
fields generated by a current in the trapping wire and 
a homogeneous external bias field i?bias- An ensemble 
of TV ^ 1.5 X 10® atoms at Tatom ~ IpK is loaded from 
the optical tweezers into this superconducting microtrap, 
formed at 400 pm from the chip surface. After adiabatic 
compression, the cloud is evaporatively cooled to achieve 
either a thermal cloud or a nearly pure Bose-Einstein con¬ 
densate (BEC). The ensemble is then magnetically trans¬ 
ported to a position z ^ 18 pm below the superconduct¬ 
ing ring by rotating Rbias around the x axis (FigQb)) 
and adjusting the current in the wire. The longitudinal 
position of the cloud along x is controlled by an addi¬ 
tional field i?conf created by a confinement wire on the 
backside of the chip, see [S] for details. 

The macroscopic superconducting ring exhibits quan¬ 
tum behavior that impacts the cold atomic cloud. In the 
superconducting state the fluxoid is quantized [35] as 

n • 4>o = PoAl^j- ds-f 4>. (1) 

This follows from the fact that the single-valuedness of 
the wavefunction requires any closed integral over the 
wave vector to be an integer multiple of 27r. Here, 
$0 = h/2e is the magnetic flux quantum and the right 
hand side needs to be evaluated along a closed contour 
within the superconductor. Al (~ 100 nm for our Nb thin 
films) is the London penetration depth, j is the super¬ 
current density and $ is the total magnetic flux through 
the closed contour. If the superconductor is large com¬ 
pared to Al, which is the case for our geometry at tem¬ 
peratures well below the transition temperature the 
integral over j can be neglected. Then, $ is quantized in 
multiples of 

n • $0 = $ = y 5 • dA (2) 

$ is given by the sum of the flux applied above Tc, 
‘hfreeze = / -Bfreeze ’ dA, and the flux LJ Created by the 
supercurrents J circulating around the ring 

7Z • 4*0 ~ 4*freeze LJ. (3) 

Here, L is the inductance of the ring and Hfreeze is the 
magnetic field applied to the ring during cooling. After 
cooling through Tc the value of n is defined as the integer 
closest to ^freeze/^o- 

Using $0 = ABfj-egze'^riro [SS], we expect for our ge¬ 
ometry a field increment ABfreeze of about 66.5 mG to 
change the flux in the ring by 14>o. Having turned off 
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FIG. 2. (a) Cross section along the dashed line in Fig. 1(a) 
with the principal magnetic x field components of the trap and 
the ring in longitudinal direction, (b) Isopotential plot of the 
calculated trapping potential for 4 flux quanta in the ring. 
An asymmetric potential with two local minima is created, 
with the lower minimum (dimple) above the ring structure 
(black markers on x axis). Each contour line corresponds to 
an energy change of /sb-SOuK. (c) Calculated potential along 
the longitudinal axis (black line Fig.j^b)) with 4 flux quanta 
in the ring (solid red) and the unperturbed harmonic trap 
(dashed). 


^freeze, the (quantized) flux through the ring is con¬ 
served by the induced circulating current Jfreeze- Any 
fields applied to the ring in the superconducting state, 
for instance by the magnetic trap, are compensated by 
screening currents Jscreen, so that the total current is 

J Jfreeze “b Jscreen • 

The magnetic (dipole) field Sring created by currents J 
locally modifies the magnetic trapping potential for the 
atoms in the vicinity of the structure. To estimate the 
contribution of Bdng and its impact on the trapping po¬ 
tential, let us consider a cigar-shaped harmonic trap with 
oscillation frequencies uix ^ ^y,z, whose radial axis y is 
centered above the ring and whose size is on the order 
of the ring diameter. The offset field B^ at the mini¬ 
mum of the trap is considered to point along the x-axis. 
The x-component of H^ing increases B^ on one side of 
the ring and reduces it on the other (Fig. [^a)), which 
leads to an asymmetric double well potential for the cold 
atomic cloud (Fig. [^b) and (c)). Hence, Bring leads to a 
position shift of the minimum of the magnetic trapping 
potential along the longitudinal axis. In addition, the 
position of the potential minimum is shifted towards the 
surface with increasing number of flux quanta in the ring 
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FIG. 3. (a) Measured integrated density profile for different 
freezing fields. Each column represents the integrated den¬ 
sity profile averaged over 9 absorption images. Adjacent lines 
differ by 7 mG in the field applied during cooling. The split 
in the density profile shows the emergence of the double well 
potential due to the ring field, (b) Relative atom number of 
an ensemble trapped at the superconducting ring, obtained 
by integrating the density profile along the a:-axis shown in 
(a). The mean atom number is calculated from 9 pictures 
per frozen field applied during cooling. The dashed vertical 
lines have a spacing of 65.9 mG, which is the measured value 
for one flux quantum. The atom number is normalized to 
the maximum number measured in the trap. The red dashed 
lines indicate the calculated atom numbers obtained from the 
numerical simulation of the Boltzmann distribution, (c) Inte¬ 
grated density profile calculated from the predicted trapping 
potential using Boltzmann distributed atoms. 


and the barrier height (trap depth) between the dimple 
and the surface is reduced. This effect leads to a decrease 
in the number and temperature of atoms trappable in this 
dimple. Using the method described in m for the simu¬ 
lation of the supercurrent densities in the trapping wire 
and the ring, we numerically calculate the held distribu¬ 
tion composed of the ring field and the trapping held with 
Biot-Savart’s law. The results of these calculations are 
in good agreement with simulations based on inductance 
calculations (3D-MLSI software package) [35]. Below the 


ring structure, the ring held for l$o leads to a held shift 
of ~ 3 mG as compared with the unperturbed harmonic 
trap (Fig.j^c)), corresponding to a dimple with a depth 
on the order of 100 nK. Furthermore, the alteration of 
the potential landscape caused by the circular supercur¬ 
rents leads to a longitudinal center-of-mass oscillation fre¬ 
quency that depends on the number of hux quanta n in 
the ring. 

For the measurements, we hrst prepare the hux state 
of the ring by heating up the chip to a temperature above 
Tc and subsequently cool it to T = 4.2 K in a homoge¬ 
neous magnetic held Breeze, applied perpendicular to the 
surface and ranging from —500 mG to 500 mG. The mag¬ 
netic helds are calibrated by microwave spectroscopy, i.e., 
the atoms are prepared in the state F = l^mp = —1 and 
the number of atoms in the state F = 2^ mp = 0 is mea¬ 
sured after application of a microwave pulse with variable 
frequency. Limited by fluctuations of the magnetic held 
in the laboratory, the absolute value of i?freeze is known 
within ±5 mG. 

A thermal cloud of Af ~ 2 x 10® atoms with a temper¬ 
ature of ^ 250 nK is prepared below the trapping wire 
and brought close to the ring, where it is held for 1 s. 
For each value of Breeze we take nine absorption images 
in situ by reflection imaging |33| along the y direction. 
After averaging over the images, we integrate the cal¬ 
culated column density along the z direction to obtain 
a one dimensional proflle of the atomic cloud along the 
axis of weak conhnement to reveal the impact of the ring 
field along x. In Fig. [^a) the density profiles are plotted 
vs. Bfreeze- We observe steps in the integrated density 
profile occurring when the number of flux quanta in the 
ring changes. For certain flux states, two distinct den¬ 
sity peaks, which indicate the double well potential, are 
discernible. By further integration of the profiles shown 
in Fig. [^a), we obtain the atom number N as a, func¬ 
tion of Bfreeze- The atom number is normalized to the 
maximum number measured in the trap and plotted in 
Fig-ib). There are clearly visible equidistant steps with 
a width of AB^eeze = 65.9 ± 2.3 mG, indicated by the 
blue vertical lines. The theoretically predicted value of 
ABfreeze = 66.5 mG per flux quantum is well within the 
error bars of the measurement. In Fig. |^b) it is visi¬ 
ble that we achieve a resolution better than single flux 
quanta. 

The measurements in Fig. Hb) are not symmetric 
around the value Bf^eeze = 0. As the magnetic trap itself 
has a magnetic field component perpendicular to the sur¬ 
face, a screening current Jscreen is induced in the ring to 
compensate this field. The screening current contributes 
to the trapping potential even for Breeze = 0. Only if the 
sum of the fields perpendicular to the ring is equal to the 
number of trapped flux quanta, the harmonic trap is un¬ 
perturbed by the screening currents. In this case, there is 
no net current around the ring, i.e. Jscreen = — Jfreeze, and 
only Meissner currents, which keep the superconducting 
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FIG. 4. Trapping frequencies measured for different freezing 
fields. The dashed vertical lines indicate the jump in the flux 
quanta. The black dots with errorbars were obtained from 
the measurement. The dotted horizontal lines are calculated 
values for different numbers of frozen flux quanta. 


film itself field free, are present m- 

To gain a qualitative understanding of the impact of 
the screening currents on the density profile, we set up 
a simplified model and numerically calculate the modifi¬ 
cation of the trapping potential for different numbers of 
flux quanta in the ring. We simulate the atomic den¬ 
sity in the trap using a Maxwell-Boltzmann distribu¬ 
tion for the energy, assuming a maximum temperature 
of Tmax = 230 nK. As the barrier height between the 
trap and the surface (trap depth) depends on the num¬ 
ber of flux quanta n in the ring, we truncate the en¬ 
ergy in the Boltzmann distribution to the trap depth, 
which changes the volume occupied by the cloud [401142] . 
Furthermore, we incorporate the loss of atoms accord¬ 
ing to a heuristic scaling of the atomic density, assuming 
p(n)/pmax = Tdepth(?T-)/Tdepth(0). The simulations take 
into account the existence of a double well potential, in 
which both traps have different trap depths. To compare 
the calculations with the observed density profiles, we 
sum over the calculated atomic density along the y and 
z directions and plot the result vs. the number of frozen 
flux quanta n. The result is shown in Fig. p c), which 
closely resembles the experimental data in Fig. ^a). The 
number of trapped atoms is estimated by additionally 
summing up the density distribution (Fig. i c)) along 
the x-direction, which leads to the red dashed lines in 
Fig. §h). Our simple model qualitatively and moder¬ 
ately quantitatively matches the behavior of the mea¬ 
sured atom number for different flux quanta n and jus¬ 
tifies the assumed scaling in the densities. We attribute 
the discrepancies between experiment and simulation to 
losses during the loading process of the double well, which 
are not taken into account in the calculations. 

In order to obtain additional information on the num¬ 
ber of frozen flux quanta, we have performed a measure¬ 
ment of the center-of-mass oscillation frequency of the 
trap at the ring for various values of Breeze • For this 
measurement, we prepared a BEG filling only the low- 


lying dimple of the potential with atoms, where the trap 
frequency is expected to vary significantly with the num¬ 
ber of flux quanta. To measure the frequency, a center- 
of-mass oscillation of the atoms along the longitudinal 
axis of the trap was excited by rapidly displacing the 
minimum of the magnetic potential using the current in 
the confinement wire. After a variable hold time (0 to 
200 ms), a microwave pulse was applied to transfer the 
atoms into the untrapped F = 2^mp = 0 state and the 
oscillation frequency was extracted from the position of 
the cloud after a time-of-fiight of 12 ms. Fig. [^ shows the 
measured oscillation frequencies along the expected val¬ 
ues extracted from our simulation of the potential. The 
dotted vertical lines are based on an atom number mea¬ 
surement similar as in Fig. [^b) and show the expected 
values of iJfreeze at which the number of flux quanta in the 
ring changes. The simulations are in qualitative agree¬ 
ment with the measurement and show our resolution on 
the level of single flux quanta; deviations between exper¬ 
iment and calculation can be attributed to uncertainties 
in the applied magnetic fields. 

In summary, we have demonstrated that a cold atomic 
cloud of ®^Rb atoms positioned close to a superconduct¬ 
ing ring is sensitive to the magnetic field created by sin¬ 
gle flux quanta. The modification of the trapping po¬ 
tential by this field is detectable in two trap characteris¬ 
tics, firstly in the trap depth and therefore in the atom 
number of the ensemble, and secondly in the trapping 
frequency inside the created dimple trap. We also ex¬ 
pect that the variation of the number of flux quanta in 
the ring will impact internal degrees of freedom of the 
trapped atoms, such as the energy difference between 
Zeeman sublevels, which is accessible by means of Ram¬ 
sey interferometry. This sensitivity paves the way to¬ 
wards future experiments that interface, e.g., SQUIDs 
and cold atomic clouds and exploit the atomic ensem¬ 
ble as a robust quantum memory |43|. Strong coupling, 
as demonstrated for example with Nitrogen vacancy cen¬ 
ters |44j, could be achieved for cold atoms by reducing 
the loop size, resulting in larger fields per flux quantum, 
or by bosonic enhancement of the coupling strength be¬ 
tween two macroscopically populated atomic states [33] . 
The impact of the applied magnetic fields, trapped flux 
within the superconducting structures, and stray light on 
the coherence of the superconducting circuits is yet to be 
studied, but is not expected to fundamentally limit the 
coupling between the two systems. 
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